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ABSTRACT 

The detection of point sources in Cosmic Microwave Background maps is usually 
based on a single-frequency approach, whereby maps at each frequency are filtered 
separately and the spectral information on the sources is derived combining the results 
at the different frequencies. On the contrary, in the case of multi-frequency detection 
methods, source detection and spectral information are tightly interconnected in order 
to increase the source detection efficiency. 

In this work we apply the matched multifilter method to the detection of point 
sources in the WMAP 7yr data at 61 and 94 GHz. This linear filtering technique takes 
into account the spatial and the cross-power spectrum information at the same time 
using the spectral behaviour of the sources without making any a priori assumption 
about it. We follow a two-step approach. First, we do a blind detection of the sources 
over the whole sky. Second, we do a refined local analysis at their positions to improve 
the signal-to-noise ratio of the detections. At 94 GHz we detect 129 5c objects at 
\b\ > 5° (excluding the Large Magellanic Cloud region); 119 of them are reliable 
extragalactic sources and 104 of these 119 lie outside the WMAP Point Source Catalog 
mask. Nine of the total 129 detections are known Galactic sources or lie in regions of 
intense Galactic emission and one additional (weak) high-Galactic latitude source has 
no counterpart in low-frequency radio catalogues. Our results constitute a substantial 
improvement over the NEWPS-3year catalogue. 

Key words: methods: data analysis - techniques: image processing - radio contin- 
uum: galaxies - cosmic microwave background - surveys 



1 INTRODUCTION 

The study of the Cosmic Microwave Background (CMB) 
provides a very useful tool to understand the Universe and 
its evolution. A proper analysis of this primordial radia- 
tion allows us to discriminate between different evolution- 
ary models of the Universe. Different experiments have taken 
data with different observing conditions (i.e. resolutions, fre- 
quencies, fields of view, ...) in order to improve our compre- 
hension of this radiation and, therefore, of the Universe. The 
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excellent sensitivities of modern instruments, close to funda- 
mental limits, imply that our ability to accurately measure 
the CMB temperature anisotropics is limited by the con- 
tamination of CMB maps by astrophysical emissions (fore- 
grounds). For this reason, many techniques have been devel- 
oped in order to separate the different components that one 
can find in the maps. While on moderate to large angular 
scales the main contaminants are diffuse Galactic emissions, 



both in temperature ( Toffolatti et all 19981: iDe Zotti et al.l 


1 19991: iHobson et all 


19991: IDe Zotti et alj|2005h and in polar- 


isation dTucci et al 


20041. 20051: Looez-Canieeo et aljboOSl; 


lAreheso et alj|201l£ 


)• 
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Astrophysical foregrounds can be both a disturbance 
for CMB studies and interesting per se. The WMAP 
surveys have made possible to investigate for the first 
time the statistical properties of bright radio sources 
above 



20031; iHinshaw et all 


2007|; |L6pez-Caniego et al.l 


2007; 


Chen & Wrightl 200S 


; iGonzalez-Nuevo et al. 


2008; 


Massardi et all 120091; 


IWright et al. 20091; Gold et al. 


201lf). The first Planck results (Planck Collaboration et al. 



201181 ) have already offered the possibility of extending 



the study to higher frequencies a nd fainter flux densities 
(|Planck Collaboration et alJl20lTbh . 

There is a rich literature on me thods to extract 
comp act sources from CMB maps (see iHerranz fc Vielval 
l20ld . for a recent review). Success ful techniqu es ex- 
ploit a variety of tools: w a velets 
2003; IGonzalez-Nuevo et all l200fj ; 



(IVielva et al.l 
janz et alj 



Lopez-Caniego et al.ll2007f). matched f ilters ( Vikhlinin et al 



dSanz et all [20oH; IChiaiw et 
Herranz et alj l2002aj; lLopez-Caniego et alj 



filtering 



1995 ; Tegmark fc de Oliveira-Costal Il998l ; Barreiro et al 
2003 ; lLopez-Caniego et al. l2006h and other linear 
tools 



Lopez-Caniego et al 



al; ILopes 
T l2005bl ). 



2004, 



12002 
12005a 



Usually these methods filter 
maps at each frequency separately. In addition, Bayesian 
techniques that include prior information about the 
distribution of the sources have been p roposed in the 
literature dHobson fc McLachlanl 120031 ; ICarvalho et al.l 
120091 ; I Argiieso et al.ll2011ri ). 

Although single-frequency filtering techniques have 
been remarkably successful, a multi-frequency approach al- 
lows us to take advantage of additional information such 
as the cross-power spectrum of the noise. This makes pos- 
sible to improve the significance of sources seen at differ- 
ent frequencies but with relatively low signal-to-noise ratio 
in each frequency channel without any a priori assumption 
on their spectral properties. However, mult i- frequency de- 
tection of p oint sources in CMB m aps is still a poorly ex- 
plored field. IHerranz fc Sanj (|2008l ) introduced the matrix 
filters technique (MTXF) as the first fully multi-frequency, 
non-parametric, linear filtering technique that is able to find 
point sources and to do unbiased estimations of their flux 
densities. The MTXF technique exploits the distinctive spa- 
tial behaviour of these s ources without assum ing any spe- 
cific spectral behaviour. IHerranz et al.1 |2009T l applied the 
MTXF to realistic simulations of the Planck radio channels, 
showing that it is possible to practically double the number 
of detections for some of the channels with respect to the 
single-frequency matched filter approach for a fixed reliabil- 
ity level. 

The MTXF approach exploits the multi-frequency in- 
formation only on the diffuse components (like CMB, Galac- 
tic emissions and noise). A more complete approach should 
take into account also the correlations among the differ- 
ent frequency maps due to point sources thems e lves. A 
step in this direction was done by IHerranz" et al.1 (|2002bl ) 
who presented the matched multi-filter (MMF) method, 
a generalisation of the standard matched filter for multi- 
frequency data where the frequency dependence of the signal 
is known. T his is the case for the thermal Suny aev-Zerdovich 
(SZ) effect (|Sunvaev fc Zeldovichlll970l . ll972r i. for which the 
method was originally developed and that has been used to 
build the Planck Early SZ catalog. Unfortunately this is not 



the ca se for extragalactic radio sources. However. lLanz et al.l 
(|2010T ) took advantage of an analogy between the SZ and the 
point source case. In the SZ case, the spectral behaviour of 
the sources is known but the size of the clusters is not. To 
deal with this, the size of the source is described, in the de- 
sign of an MMF, by means of a free parameter, whose value 
is determined m aximising the signal-t o-noise ratio for each 
detected source (|Herranz et al.ll2002bl ). In the case of point 
sources, the angular profiles are known (they are given by 
the point spread functions, PSFs, of the instrument) and 
we only need to optimise a set of parameters describing the 

spectra l shape. 

In lLanz et all (|2010T ) the MMF was applied to realis- 
tic simulations of the Planck mission. The promising results 
obtained in that work prompted us to apply the method 
to real data, namely to the 7-yr WMAP data. More pre- 
cisely, we have considered the V and W maps (61 and 94 
GHz respectively), because at these frequencies the knowl- 
edge of the statistical properties of radio galaxies was poor. 
The structure of the article is as follows. In §[2] we briefly 
describe the method, and in §[3] the data used. In § [J] we 
present and discuss our results, comparing our flux density 
estimates with measurements with other instruments. Our 
main conclusions are summarised in §[S] 



2 METHOD 

The MMF is the optimal linear detection method when the 
spatial profile and the frequency dependence of the sources 
are known. By 'optimal' we mean, as it is common in statis- 
tics, that the estimation of the flux density of the sources is 
unbiased and has minimum variance (maximum efficiency). 
In Fourier space the MMF takes the form: 



*(q) = a P _1 F, 



1 = / dq F T P _1 F, 



(1) 



where q is the Fourier mode, ^f(q) is the column vector of 
the filters *&(q) = [tp^(q)], F is the column vector F = [f v T v ], 
f v is the frequency dependence, T v (q) is the source profile at 
each frequency v and P _1 is the inverse matrix of the cross- 
power spectrum P. The MMF takes as arguments N images 
(in this work N = 2) and returns a single filtered image 
where the sources are optimally enhanced with respect to 
the noise. The variance of the output filtered image is given 
by: 



dq<J> P* = a. 



m 



If a source has a signal-to-noise ratio ^ 5 after the filtering, 
we say that we have a detection at that position. 

Since we are dealing with point sources, the profiles t„ 
are directly given by the PSFs of the instrument. The cross- 
power spectrum P is not known a priori but can be inferred 
from the data under the assumption that the point sources 
are sparse. The unknown frequency dependence, f = [fv], of 
a given source can be parametrised as: 

Sy _ ( v 
So \ Vq 

where S v and So are the flux densities at a frequency v and 
at the frequency of reference v p . The spectral ind ex 7 is a free 
parameter. As pointed out bv lLanz et al.1 (|2010T l. each image 



(3) 
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is filtered several times for different MMFs. These MMFs are 
identical except for the spectral index 7. The test values of 
the spectral index used in this work are —3.5 ^ 7 ^ 3.5, with 
a step of 0.05. As it was shown in the cited work, the signal- 
to-noise (SNR) ratio of the detected source is maximal for 
the correct choice of this parameter, and by construction (we 
want an efficient estimator of the flux density of the source), 
the uncertainty assigned to the source is the square root of 
the variance expressed in eq. @ for the correct value of 7. 
While the parametrisation of eq. ((3} is perfectly adequate in 
our case, any other parametrisation, even a non-functional 
description of the vector f by means of its components, could 
be used in other cases, e.g. when more than two frequencies 
are considered simultaneously. 



3 APPLICATION TO WMAP DATA 

As mentioned above, the WMAP V and W bands are par- 
ticularly interesting because only a small fraction of WMAP 
sources were detected with SNR ^ 5a in these bands. We 
have used an adapted ver s ion of the code desc ribed in 
120071) . iMassardi et"al] (|2009l ) and in 
2003 ) , modified in order to handle two 
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frequencies simultaneously and to accommodate our MMF 
filtering. The code reads in an input parameter file contain- 
ing the specific characteristics of the maps to be analysed 
as well as the patch size, the pixel size and overlap among 
the patches to effectively cover the 100% of the sky. Then, 
it reads in the two input maps in FITS format and ex- 
tracts the patches to be analysed using the tangential plane 
approximation implemented in the CPAC library [_. Each 
pair of V and W patches is analysed simultaneously using 
the MMF, and a first set of detections is produced. Next, 
the code iteratively explores different values of the spec- 
tral index for each source, allowing for the appearance of 
new possible detections. At the end, the code gathers all 
the candidates and generates a list of detections above a 
given SNR, converting the positions of the detected objects 
in the plane to the sphere. We have used flat patches of 
14.6° x 14.6°, each containing 128 x 128 pixels. The pixel 
area is 6.871' x 6.871', correspondi ng to the HEALPix res- 
olution parameter N ai d e = 512 fsee lGorski et all (|2005h for 
more details). In total we have 371 patches, with significant 
overlaps among them in order to a void as much a s poss i- 
ble border effects. As described in iMassardT et al.l (|2009h . 
we perform a two-step process, first doing a blind search 
across the sky and then refining the analysis obtaining new 
patches centred at the positions of the sources identified in 
the blind step. In this way we further reduce any projection 
and border effect that we could have. 

When dealing with the V and especially the W WMAP 
bands it is necessary to carefully characterise the properties 
of the beams, both to build the optimal filters (eq. _) and 
to obtain a good photometry of the sources. WMAP beams, 
particularly at 94 GHz, are highly non-Gaussian and non 
circularly symmetric. We also need to take into account that 
for the W band the pixel size (6.871') is dangerously close 
to the minimum sampling of the FWHM required to avoid 



aliasing, and that real sources are not necessarily located 
at the geometrical centre of the pixels. We have dealt with 
these systematic effects in the following way. 

We used the beam transfer functions (circularly sym- 
metric) provided by WMAP_| to create high resolution tem- 
plates of the beams (equivalent to the HEALPix N si d e = 
4096) projected into the tangent plane. Please note that the 
maps that we have used to perform our analysis have a pixel 
size corresponding to N s id e = 512. In addition, we must 
consider that sources are not necessarily located at the geo- 
metrical centre of these pixels. Therefore, taking the beam 
transfer functions, we divided each N si d e = 512 pixel in 8 x 8 
subpixels (N a id e = 4096) using the following formula: 



fe; =2tt / b s {6)Pi(cos6)d(cos6)/Q,L 



(4) 



where we obtained the value of the beam template in any an- 
gular position (6) from the beam transfer function using the 
Legendre polynomials. Then we averaged the high-resolution 
beams over all the 8x8 possible displacements of the source 
centre inside a N 3 id e = 512 pixel. In order to obtain a well- 
estimated average of the beams, we repeated this process 
5000 times. These averaged beam templates were then de- 
graded to the same resolution of our maps (N a id e = 512). In 
this way we take into account, on average, the effects of pix- 
elisation and the possible offsets of the sources with respect 
to the geometrical centre of the pixel. 

A careful correction for the non-circularity of the beams 
would require a detailed knowledge of the combination 
of beam orientations for all scans through every posi- 
tion in the sky. Since this information is not available to 
us, we only applied to the flux densities an average cor- 
rection factor that takes into account the effective beam 
areas for each channel. The correction factors were ob- 
tained using the self-calibration method based on SEx- 
tractor jBertin fc Arnouts! 1 19961) AUT O photometry, de- 
scribed in iGonzalez-Nuevo et al.l (|2008l ). In that work, the 
NEWPS 3yr sources detected above the 5a level by SEx- 
tractor and with flux densities S ^ 1 Jy were used for the 
self-calibration. Only 3 W-band sources satisfied these re- 
quirements and therefore the correction factor for this band 
could not be reliably estimated. To overcome this problem 
we have applied SExtractor to the W map at the positions 
of 296 sources observed by the Australia Telescope Com- 
pact Array (ATCA) and the Very Large Array (VLA). Fif- 
teen of them were detected at ^ 3a by SExtractor. Since 
all the ATCA/VLA sources have been detected at a very 
high significance level, we were confident that there were 
no spurious detections in this sample. These sources were 
used to recalculate the effective beam area at 94 GHz. We 
obtained (2.50 ± 0.09) x 10~ 5 sr, to be compared with the 
nominal value of 2.07 x 10~ 5 sr for the symmetrised beam. 
Note that all the ingredients described in this section must 
be introduced in the code before the filtering process. 
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Figure 1. Position on the sky, in Galactic coordinates, of the 129 sources detected by the MMF on the at 61 and 94 GHz WMAP 7yr 
maps. Blue crosses represent the position of the 119 sources with confirmed counterparts in other catalogues (see Table |2j; red squares 
represent the positions of the 10 Galactic or unidentified detections (see Table [TJ. A coloured version of this figure is available with the 
online edition of the paper. 
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RA (°) Declination (°) S94GHZ (Jy/beam) rms (Jy) spectral index Id in other catalogues 



13.237 


56.576 


2.86 


0.22 


-1.25 


NGC 281 


44.622 


-18.834 


0.54 


0.10 


-2.40 




52.264 


31.338 


4.81 


0.57 


4.50 


NGC 1333 


60.731 


36.189 


0.49 


0.09 


-3.35 


NGC 1499 


83.654 


22.045 


155 


2 


-1.15 


Crab Neb 


83.807 


-5.419 


216 


3 


-0.35 


Orion region 


84.018 


-6.338 








Orion region 


85.396 


-1.917 


30 


2 


-0.2 


Orion region 


246.756 


-24.633 








Ophiucus region 


285.398 


-36.983 


3.31 


0.55 


3.70 


NGC 6729 



Table 1. Equatorial coordinates and identifications of the Galactic or unidentified detections. Sources without data are detections that, 
due to the multifrequency method, do not have a clear estimation of their flux density of the spectral index inside the studied range of 7. 
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Table 2: Sources detected by the MMF on WMAP 7yr maps with coun- 
terparts in lower frequency catalogues. Presumably Galactic sources are 
listed in Table [1] The Equatorial coordinates (first two columns) are in 
degrees. The flags in the last column indicate if the source is listed in 
any o f the the WMAP 7yr catalogues at any frequency (W) ( Gold et al] 
120111 ) . in NEWPS-5yr at 61 GHz (N) jMassardi et~al]|200g) or in the 
Planck catalogue at 70 or 100 GHz (P) ( Planck Collaboration et alj 
l2011aT ). Sources without data are detections that, due to the multifre- 
quency method, do not have a clear estimation of their flux density of 
the spectral index inside the studied range of 7. These sources with no 
flux density information are not taken into account in the subsequent 
analysis. 



RA O 


Declination (°) 


Smghz (Jy/beam) 


rms (Jy) 


spectral index 


Id in other catalogues 


Flags 


1.483 


-6.350 


1.30 


0.23 


-0.80 


PMN J0006-0623 


WNP 


5.033 


73.493 


1.41 


0.26 


-0.45 


GB6 J0019+7327 


P 


6.514 


-35.154 


0.72 


0.14 


-1.65 


PMN J0026-3512 


WN 


9.603 


-24.892 


1.65 


0.33 


0.30 


PMN J0038-2459 


WNP 


16.715 


-40.600 


1.74 


0.21 


-0.55 


PMN J0106-4034 


WNP 


23.160 


-16.928 


0.91 


0.17 


-1.20 


PKS 0130-17 


WNP 


24.279 


47.846 


1.91 


0.19 


-1.20 


GB6 J0136+4751 


WNP 


24.383 


-24.518 


1.56 


0.30 


0.20 


PMN J0137-2430 


WNP 


32.749 


-51.070 


1.93 


0.17 


-1.00 


PMN J0210-5101 


WNP 


34.293 


73.803 


0.73 


0.13 


-1.70 


GB6 J0217+7349 


NP 


39.510 


28.836 


2.20 


0.33 


-0.20 


GB6 J0237+2848 


WNP 


39.723 


16.654 








GB6 J0238+1637 


WP 


43.374 


-54.666 


1.56 


0.20 


-0.55 


PMN J0253-5441 


WNP 


46.088 


-62.233 


1.31 


0.25 


0.05 


PMN J0303-6211 


WNP 


49.957 


41.450 


2.67 


0.16 


-1.90 


GB6 J0319+4130 


WNP 


50.377 


-37.092 


0.81 


0.14 


-1.55 


Fornax A 


NP 


53.549 


-40.203 


1.47 


0.25 


-0.20 


PMN J0334-4008 


WNP 


54.252 


32.327 


2.57 


0.42 


0.70 


GB6 J0336+3218 


P 


54.862 


-1.734 


2.33 


0.33 


0.05 


PMN J0339-0146 


WNP 


57.187 


-27.832 


1.36 


0.26 


0.00 


PMN J0348-2749 


WNP 


61.062 


-36.120 


2.14 


0.21 


-0.65 


PMN J0403-3605 


WNP 


64.548 


38.025 


1.96 


0.16 


-2.15 


3C 111 


NP 


65.845 


-1.374 


3.39 


0.25 


-1.15 


PMN J0423-0120 


WNP 


67.166 


-37.869 


1.66 


0.29 


0.05 


PMN J0428-3756 


WNP 


68.301 


5.418 


3.51 


0.41 


0.50 


GB6 J0433+0521 


WNP 


72.260 


-81.053 


1.61 


0.28 


0.25 


PMN J0450-8100 


WNP 


73.914 


-46.248 


3.05 


0.26 


-0.50 


PMN J0455-4616 


WNP 


74.270 


-23.427 


1.66 


0.20 


-0.75 


PMN J0457-2324 


WNP 


78.849 


-45.994 


1.02 


0.19 


-1.40 


PMN J0515-4556 


WP 


79.975 


-45.795 


2.30 


0.21 


-1.20 


PMN J0519-4546 


WNP 


80.849 


-36.456 


2.57 


0.19 


-1.05 


PMN J0522-3628 


WNP 


84.678 


-44.028 


3.15 


0.21 


-1.15 


PMN J0538-4405 


WNP 


91.968 


-6.396 


7.99 


0.49 


0.40 


PKS 0605-06 


NP 


92.460 


-15.756 


1.79 


0.24 


-0.50 


PMN J0609-1542 


WNP 


97.363 


-19.946 


1.66 


0.32 


0.55 


PMN J0629-1959 


WNP 


98.985 


-75.242 


2.03 


0.19 


-0.80 


PMN J0635-7516 


WNP 


101.703 


44.896 


1.80 


0.33 


0.05 


GB6 J0646+4451 


WNP 


102.564 


-16.699 


1.75 


0.26 


-0.35 


PKS 0648-16 


NP 


110.556 


71.313 


2.15 


0.26 


-0.05 


GB6 J0721+7120 


WNP 


111.493 


-0.972 


2.23 


0.38 


0.65 


PMN J0725-0054 


WP 


114.856 


1.724 


0.96 


0.18 


-1.35 


GB6 J0739+0136 


WNP 


117.689 


12.478 


2.04 


0.23 


-0.90 


GB6 J0750+1231 


WNP 


126.524 


3.272 


1.27 


0.22 


-0.90 


GB6 J0825+0309 


WNP 


129.076 


-20.280 


1.07 


0.19 


-1.10 


PMN J0836-2017 


WNP 


130.356 


70.885 


0.91 


0.15 


-1.30 


GB6 J0841+7053 


WNP 


133.730 


20.128 


2.48 


0.29 


-0.70 


GB6 J0854+2006 


WNP 


140.233 


44.738 


0.89 


0.14 


-2.00 


GB6 J0920+4441 


WP 
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Table 2 — continued 



RA (°) 


Declination (°) 


594GHz (Jy/beam) 


rms (Jy) 


spectral index 


Id in other catalogues 


Flags 


140.441 


-26.324 


1.07 


0.21 


-0.85 


PMN J0921-2618 


WP 


141.810 


38.974 


2.63 


0.21 


-1.15 


GB6 J0927+3902 


WNP 


159.372 


-29.584 


1.29 


0.22 


-0.75 


PMN J1037-2934 


WNP 


159.659 


5.136 


0.80 


0.16 


-1.70 


GB6 J1038+0512 


W 


164.323 


-80.083 


1.65 


0.19 


-0.85 


PMN J1058-8003 


WNP 


164.615 


1.562 


2.86 


0.22 


-1.10 


GB6 J1058+0133 


WNP 


171.870 


-18.836 


0.97 


0.19 


-1.20 


PMN J1127-1857 


WP 


176.717 


-38.190 


0.92 


0.15 


-1.75 


PMN J1147-3812 


WNP 


178.285 


49.483 


1.47 


0.19 


-0.70 


GB6 J1153+4931 


WNP 


179.831 


29.192 


1.21 


0.22 


-0.60 


GB6 J1159+2914 


WNP 


187.270 


2.104 


10.71 


0.33 


-1.20 


GB6 J1229+0202 


WNP 


187.699 


12.338 


5.99 


0.26 


-1.05 


GB6 J1230+1223 


WNP 


191.732 


-25.848 








PKS 1244-255 


WNP 


193.992 


-5.821 


8.59 


0.32 


-1.15 


PMN J1256-0547 


WNP 


194.823 


51.645 


1.37 


0.25 


0.05 


GB6 J1259+5141 


W 


197.630 


32.472 


0.82 


0.13 


-1.75 


GB6 J1310+3220 


WNP 


201.390 


-42.960 


16.92 


0.46 


-1.20 


Centaurus A 


NP 


204.426 


-13.021 


4.96 


0.29 


-0.60 


PMN J1337-1257 


WNP 


209.267 


19.351 


1.39 


0.21 


-0.75 


GB6 J1357+1919 


WNP 


214.876 


54.421 


1.05 


0.20 


-0.50 


GB6 J1419+5423 


WNP 


216.967 


-42.070 


1.75 


0.23 


-0.75 


PMN J 1427-4206 


WNP 


229.382 


-24.370 


2.31 


0.28 


-0.40 


PMN J1517-2422 


WNP 


237.415 


2.563 


1.54 


0.27 


-0.55 


GB6 J1549+0237 


WNP 


237.763 


5.461 


1.77 


0.26 


-0.50 


GB6 J1550+0527 


WNP 


243.420 


34.145 


2.58 


0.25 


-0.30 


GB6 J1613+3412 


WNP 


244.386 


-77.279 


1.07 


0.16 


-1.20 


PMN J1617-7717 


WNP 


246.561 


-29.832 


1.91 


0.30 


-0.40 


PKS 1622-29 


NP 


248.816 


38.155 


3.25 


0.25 


-0.50 


GB6 J1635+3808 


WNP 


249.747 


57.392 


1.98 


0.29 


0.55 


GB6 J1638+5720 


WNP 


250.551 


68.974 


0.80 


0.08 


-2.45 


GB6 J1642+6856 


WNP 


250.742 


39.825 


3.89 


0.23 


-0.75 


GB6 J1642+3948 


WNP 


260.115 


-1.025 


2.18 


0.27 


-0.40 


PKS 1717-00 


NP 


260.873 


-65.025 


1.32 


0.25 


-0.25 


PMN J1723-6500 


WP 


263.303 


-13.051 


3.05 


0.29 


-0.30 


PKS 1730-13 


NP 


263.593 


38.888 


1.08 


0.20 


-0.70 


GB6 J1734+3857 


WNP 


265.933 


-3.813 


2.55 


0.22 


-0.90 


PKS 1741-03 


NP 


267.966 


9.655 


4.05 


0.26 


-0.45 


GB6 J1751+0938 


WNP 


268.340 


28.738 


1.65 


0.21 


-0.50 


GB6 J1753+2847 


WNP 


270.301 


78.459 


1.05 


0.17 


-1.00 


NVSS J180045+782804 


WP 


270.400 


44.088 


1.40 


0.22 


-0.45 


GB6 J1801+4404 


WNP 


271.711 


69.794 


1.34 


0.20 


-0.25 


GB6 J1806+6949 


WNP 


277.443 


48.800 


1.07 


0.13 


-1.60 


GB6 J1829+4844 


WNP 


280.720 


68.162 


0.82 


0.14 


-0.90 


GB6 J1842+6809 


WP 


282.241 


67.063 


0.84 


0.13 


-1.15 


GB6 J1849+6705 


WNP 


287.737 


-20.060 


2.12 


0.34 


0.10 


PMN J1911-2006 


WNP 


290.987 


-21.081 


1.50 


0.25 


-0.65 


PMN J1923-2104 


WNP 


291.214 


-29.282 


7.01 


0.28 


-1.00 


PMN J1924-2914 


WNP 


292.018 


73.976 


1.56 


0.15 


-1.20 


GB6 J1927+7357 


WNP 


299.461 


-38.755 


2.04 


0.26 


-0.55 


PMN J1957-3845 


WNP 


299.942 


40.752 


6.02 


0.36 


-2.25 


Cygnus A 


NP 


300.337 


-17.769 


0.84 


0.16 


-1.60 


PMN J2000-1748 


WNP 


302.849 


-15.774 


1.45 


0.29 


-0.25 


PMN J2011-1546 


WNP 


314.018 


-47.226 


1.87 


0.25 


-0.45 


PMN J2056-4714 


WNP 


323.562 


-1.791 


1.25 


0.22 


-1.00 


PMN J2134-0153 


WNP 


324.128 


0.682 


1.67 


0.28 


-0.30 


GB6 J2136+0041 


WNP 


327.069 


6.999 


3.98 


0.23 


-1.15 


GB6 J2148+0657 


WNP 


327.945 


-30.455 


1.44 


0.24 


-0.55 


PMN J2151-3028 


WNP 


328.366 


47.275 


1.93 


0.26 


-0.25 


GB6 J2153+4716 


NP 


329.184 


-69.703 


1.81 


0.22 


-0.50 


PMN J2157-6941 


WNP 
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Table 2 — continued 

RA (°) Declination (°) S94GHZ (Jy/beam) rms (Jy) spectral index Id in other catalogues Flags 



oof) fi/lQ 


4Z.ZD1 


9 fifi 
Z.Oo 


u.zo 


-U.oU 


n'Rfi 79909^/191 fi 
LrDO JZZUZ^4Z10 


WTVP 
VV IN tr 


330.745 


31.724 


1.66 


0.28 


-0.05 


(jDO JzzUo+ol45 




330.851 


17.436 


1.24 


0.23 


-0.70 


GB6 J2203+1725 


WNP 


334.914 


63.283 


2.98 


0.57 


2.90 


GB6 J2219+6317 


NP 


336.405 


-5.002 


3.43 


0.26 


-0.75 


PMN J2225-0457 


WNP 


337.427 


-8.582 


2.69 


0.30 


-0.35 


PMN J2229-0832 


WNP 


338.171 


11.674 


4.19 


0.34 


0.20 


GB6 J2232+1143 


WNP 


338.909 


-48.584 


1.33 


0.19 


-0.85 


PMN J2235-4835 


WNP 


341.600 


-12.087 


1.24 


0.21 


-1.05 


PMN J2246-1206 


WP 


343.547 


16.149 


7.96 


0.35 


-0.35 


GB6 J2253+1608 


WNP 


344.514 


-28.003 


3.58 


0.28 


-0.35 


PMN J2258-2758 


WNP 


352.347 


-47.438 


1.18 


0.22 


-0.55 


PMN J2329-4730 


WP 


359.607 


-53.218 


1.24 


0.20 


-0.60 


PMN J2357-5311 


WP 
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4 RESULTS AND DISCUSSION 

The source detection was performed applying the MMF fil- 
ter to the WMAP 7yr V and W maps, at 61 and 94 GHz 
respectively. We have used only one (VI) of the 2 differenc- 
ing assemblies for the V band, because VI and V2 are very 
similar. As for the W band, we used 2 (W2 and W3) out 
of the 4 differencing assemblies because their symmetrized 
beam profiles are geometrically very similar (same as for Wl 
and W4) and this makes the photometry easier. The W2 and 
W3 maps were combined pixel by pixel, weighting with the 
inverse of the variances of the pixels. To minimise the spuri- 
ous detections due to the complex structure of the Galactic 
emissions near the equatorial plane and of the Large Magel- 
lanic Cloud (LMC) region, we have removed from the final 
catalogue objects with \b\ ^ 5° and within a radius of 5° 
around the LMC ((80.894°, -69.756°) in equatorial coordi- 
nates). No sources are detected within the Small Magellanic 
Cloud. 




Figure 2. Spectral index distribution of sources in Table [5] 



4.1 Point source detection 

The two-step filtering approach described in section [3] yields 
129 5a detections in the 61 and 94 GHz WMAP 7yr maps, 
outside the Galactic plane and LMC regions specified above. 
For each of them, the MMF gives the flux density at the ref- 
erence frequency (94 GHz, vq in eq. (|3J) and the spectral 
index, the only free parameter of the method. The uncer- 
tainty on the spectr al index can be est imated from intensive 
simulations made in lLanz et all (|201fJ ). Table[Qlists the nine 
sources that are either well known Galactic sources (like the 
Crab Nebula, the HII regions NGC 281 and NGC 1499 or 
the star-formation regions NGC 1333 and NGC 6729) or are 
located within Galactic molecular cloud complexes, like the 
Orion or the Ophiucus regions, or within regions of intense 
Galactic emissions, plus one unidentified. Table [T] gives the 
flux densities within the WMAP 94 GHz beam. In most 
cases, these are lower limits to the total flux densities be- 
cause the sources are resolved. Given the different beam sizes 
at the two frequencies, the spectral indices cannot be reli- 
ably estimated in the case of resolved sources. 

Excluding the ten sources listed in Table [T] the MMF 
has detected 119 objects, listed in Tableful All these objects 
have a counterpart in lower frequency radio catalogues, that 
cover the whole sky to much fainter flux density limit^f], 
(see Figure [T]). The apparently ultra-steep (7^ = —2.25) 
spectrum of Cygnus A is partly due to resolution effects 
(the 94 GHz beam encompasses a lower fraction of the total 
flux density than the 61 GHz beam). 

In Tables [T] and [21 sources without an estimation of the 
flux density or the spectral index are objects that are only 
'detected' at one frequency. Therefore, because the method 
we use to filter the maps is multifrequency, the spectral in- 
dices of these sources fall outside the range studied in this 
work. These sources are not taken into account in the sub- 
sequent analysis. 

The spectral index distribution of the sources in Ta- 
ble [2] shown in Figure [2] has a median value of -0.65 



3 Except for a few of the sources in Table [TJ none of the objects 
detected by the MMF has an IRAS counterpart. 



S (Jy) 


S 5 / 2 dN/dS (Jy 3/2 sr- 1 ) 




61 GHz 


94 GHz 


1.27 


10.1±1.7 


10.9±1.8 


2.01 


25±4 


18±3 


3.18 


21±5 


19±5 


5.05 


27±8 


9±5 


8.00 


14±8 


18±9 


12.68 


27±16 


9±9 


20.10 


18±18 


18±18 


31.85 


36±36 





Table 3. Euclidean-normalised differential number counts at 61 
and 94 GHz, based on our detections, in bins of A log(S) = 0.2 
and taking into account the Bayesian correction to the flux den- 
sities. 



with a dispersion of 0.71. For comparison, IWright et al.l 
(|2009l ). mostly based on lower frequency WMAP data, 
found a mean spectral index (7) = —0.09, with a 
dispersion a — 0.28. Our results are therefore con- 
sistent with the steepening of source sp ectra above 
~ 70 GHz found for different data sets dSadler et al.l 
20081; iGonzalez-Nuevo et alj I200SI; " 



Marriage et al. 



Massardi et al]|2010r . IPlanck Collaboration et al]|2011bh 



The Euclidean normalised source counts at 61 and 94 
GHz are given in Table [3] and shown in Figure [3] where they 
are c ompared with the counts in the ne arest Planck chan- 
nels (|Planck Collaboration et al.ll2011bl ) and with the pre - 
dictions of the lDe Zotti et all (120051 ) and lTucci et alj (|201lh 
models. 

In this figure, the flux densities obtained with the MMF 
were corrected with a Bayesian approach (|Herranz et al.l 
120061 ) in order to remo ve as much as possible the Eddington 
bias jEddingtonlll913h . This approach takes into account 
the distribution in flux density of the objects as a power 
law with unknown slope, and an additive Gaussian noise. 
It is important to point out that this correction is statisti- 
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61 GHz 




S(Jy) 

94 GHz 




S(Jy) 

Figure 3. WMAP 7yr Euclidean-normalised differential number 
counts at 61 and 94 GHz based on the MMF sample (red points) 
with Poissonian error bars and the Bayesian correction to the flux 
densi ties. The grey crosses are the Plan ck counts at 70 and 100 
GHz jPlanck Collaboration et alj|2011bl) . The solid lines are the 
prediction by the iTucci et al.l < 201ll) mod el at 61 and 94 GHz 
and the dashed lines the prediction by the |Pe Zotti et aL I 1120051 ) 
model at the same frequencies. 



cal, and therefore it has been taken into account only in the 
estimation of the source counts. 

The comparison in Figure [3] shows that our complete- 
ness limit is ~ 2 Jy. Above t his limit, the agree ment with the 
Planck counts and with the lTucci et alj (120111) model is gen- 
erally good. This confirms that the Tucci et alj (|201ll ) model 
deals appropriately with the high frequency behaviour of 
source spectra. 

Al though several data set s had suggested (| Sadler et al.l 

l20Qg|; iGonzalez-Nuevo et all 120081') and then detected 



( Planck Collaboration et alj 2011bl Jl a break in the bright 
extragalactic radio sources (at high-flux level) at ~ 70 GHz, 
this break is only well explaine d with the current data if we 
consider the ITucci et al.l (|201l] ) model. 



4.2 Comparison with the WMAP catalogues 

The WMAP Five-Band Point Source Catalog |GoId et al.1 
120111 ) lists all the sources that where detected at ^ 5a in 
at least one frequency channel, outside a mask excluding 
sources in the Galactic plane and Magellanic Cloud regions. 
Flux densities in the other channels are reported if the sig- 
nal is detected at a > 2a significance. In this catalogue, 94 
GHz flux densities are reported for 236 sources. However it 
is not clear how many of them are detected at a ^ 5a level 
at 94 GHz. Since the WMAP and NEWPS catalogues have 
comparable statistical completeness and reliability proper- 
ties, on the basis of the NEWPS3yr catalogue (see section 
I4.3|l we may roughly estimate that the ^ 5a detections are 
22. We expect a similar or smaller number of sources with 
the same significance in the less complete WMAP catalogue 
at the same frequency. By contrast, our method guarantees 
that all the 119 (presumably) extragalactic sources, 104 of 
which are outside the WMAP Point Source Catalog mask, 
are detected at ^ 5a at 94 GHz. 



4.3 Comparison with NEWPS catalogues 

lLopez-Caniego et al.1 |2007l ) looked for sources in the 
WMAP 3-yr maps at the positions of known 5 GHz sources 
(non-blind search) and reported 22 detections at 94 GHz 
(listed in the NEWPS3yr catalogue). iMassardi et al.1 (|2009l ) 
combined a non-blind (based on the ATCA 20 GHz survey 
catalogue) and a blind search on WMAP 5-yr maps, exclud- 
ing the W (94 GHz) channel (NEWPS5yr catalogue). At 61 
GHz they detected 169 sources, 61 of which are not present 
in our catalogue (if we also take into account Table [TJ ; 45 
of them however are detected at 3a with the MMF in 
the blind step of our detection process. Seven of the other 
18 sources are at \b\ < 15° and may be contaminated by 
Galactic emission, and three lie in the LMC region, excluded 
from our search. The remaining 8 sources have not been de- 
tected by the MMF. On the other hand, our catalogue of 
extragalactic objects contains 16 sources not present in the 
NEWPS5yr catalogue. We must remark that our catalogue 
is completely based on a blind search, whereas the NEWPS 
catalogue was guided by a previous selection of sourecs de- 
tected at low frequency. 



4.4 Comparison with ATCA and NRAO flux 
densities 

Our sample of 94 GHz detections includes 85 sources with 
ground based ~ 3 mm observa tions either with ATCA0 or 
with the NRAO 12m telescope (|Holdawav et al Il994l f1. The 
ATCA flux densities for these sources have been collected in 
the frequency range between 85 and 105 GHz with the "old" 
ATCA digital correlator with up to 2 x 256 MHz bandwidth. 
The MMF flux density estimates are compared with ATCA 
and NRAO measurements in Figure|4] There is evidence that 
the MMF somewhat overestimates the flux densities below 
~ 2 Jy, most likely due to the Eddington bias. Above 2 Jy the 



4 http:/ /www. narrabri.atnf.csiro.au/calibrators/index. html 

5 http:/ /www. alma. nrao.edu/memos/html- 
memos / almal23/memol23.pdf 
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70 GHz 




Figure 4. Comparison between 94 GHz flux densities recovered 
by the MMF and the ATCA (3 mm) and NRAO (90 GHz) flux 
density measurements, for the 85 sources in common. No correc- 
tion for the slightly different frequencies of the sets of measure- 
ments was made. The solid line is x=y. 



mean absolute fractional error (|5mmf — Sgroundl/S'ground) — 
38% somewhat higher than expected from the combination 
of nominal measurement errors and variability, suggesting 
that the true errors associated to the MMF flux density 
estimates are somewhat larger than the nominal values. In 
this respect, it must be noted that the WMAP 7-yr maps are 
averages over seven years of observations while the ATCA 
and NRAO measurements refer to a single epoch. 



4.5 Comparison with the Planck Early Release 
Compact Source Catalogue (ERCSC) 

Thanks to its higher sensitivity, Planck has detected far 
more sources than WMAP. The Early Release Compact 
Source Catalogue (ERCSC) lists 1381 sources detected at 
100 GHz and 599 detected at 70 GHz; 308 and 788 of the 
ERCSC 70 GHz and 100 GHz sources, respectively, are out- 
side of the WMAP Point Source Catalog mask. Three of 
our sources in Table [2] are not present in the ERCSC, but 
are present in lower frequency catalogues with flux densities 
consistent with those inferred from WMAP data. 

The ERCSC gives four different measures of flux den- 
sity for each source. For the comparison with our results, 
we have chosen the estimate called FLUXDET, which is 
calculated using an approach similar to ours and appears 
to have higher reliability for low signal-to-noise ratios. Fig- 
ure [5] compares the FLUXDET flux densities at 70 and 100 
GHz with the MMF ones at the nearest WMAP frequen- 
cies (61 and 94 GHz, respectively). Again there is evidence 
that the MMF flux densities are affected by the Eddington 
bias below ~ 2 Jy. The MMF flux densities of the 5 sources 
with 5mmf,94GHz ^ 8 Jy are all lower than the the ERCSC 
flux densities at 100 GHz. Two of these sources (the one at 
RA = 91.968° , dec = -6.396° and Cen A) has a strong dust 
emission, seen in the IRAS survey, that enhances the 100 
GHz flux density. The other 3 are well known highly vari- 




Figure 5. Comparison of the flux densities recovered by the MMF 
on WMAP maps at 61 and 94 GHz with the FLUXDET values 
reported in the Planck ERCSC at 70 GHz (upper panel) and 100 
GHz (lower panel), respectively. The solid line is x=y. 



able blazars (3C273, 3C279, 3C454.3), caught by Planck in 
a bright phase. 

The spectral index distributions of our sources are 
pretty close to that of the ERCSC sources: the median spec- 
tral index of our sources is -0.65, with a standard deviation 
of 0.71; for ERCSC sources with 7^2 (to avoid strong 
contamination by dust emission) and excluding the Galac- 
tic plane and the LMC region as defined previously in order 
to have a sample comparable with ours, we find a median 
value of —0.39 with a standard deviation of 0.52. 



5 CONCLUSIONS 

The detection of extragalactic point sources is a crucial task 
in the analysis of CMB maps because point sources are the 
main contaminant of the CMB power spectrum on small an- 
gular scales. From the same maps we can extract interest- 
ing astrophysical information about the point sources them- 
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selves. The development of efficient detection tools is there- 
fore very important. The MMF holds the promise of achiev- 
ing detections down to fainter flux densities than achievable 
by standard methods. 

Applying this tool to WMAP 7yr maps at 61 and 94 
GHz simultaneously we have obtained 129 5a detections at 
both frequencies, outside the Galactic plane (i.e. at \b\ > 5°) 
and the LMC regions. Nine of these sources, listed in Table[T] 
are either known Galactic sources or reside in regions of 
high Galactic emission. One additional source, also listed in 
Table [T] does not have a counterpart in low frequency radio 
catalogues. All the other 119 sources, listed in Table [2] do 
have a low frequency counterpart; 104 of them reside outside 
the WMAP Point Source Catalog mask. For comparison, 
the NEWPS-3year catalogue contains 22 5a detections at 
94 GHz. Although the Planck ERCSC contains far more 
sources at the frequencies (70 and 100 GHz) nearest to those 
used in the present analysis, we have detected three sources 
not present in the ERCSC, yet with flux densities consistent 
with lower frequency measurements. 

A comparison of the flux densities yielded by the MMF 
with those obtained at ~ 90 GHz with the ATCA or the 
NRAO 12m telescope, and with the Planck ERCSC data at 
70 and 100 GHz, shows a generally good agreement above 
~ 2 Jy, although the rms differences between MMF and 
ground based or ERCSC values are larger than expected 
from variability and nominal measurement errors. This sug- 
gests that the errors associated to MMF flux density esti- 
mates are somewhat larger than the nominal values listed 
in Table [2] Below 2 Jy, the MMF flux densities tend to be 
somewhat overestimated, as the effect of the Eddington bias. 

The distribution of 61-94 GHz spectral indices of 
sources in Table [5] has a median value eq ual to 7ef = 
—0.65 (with a dispersion of 0.71), while IWright et al.l 
(|2009l ). mostly based on lower frequency WMAP data, 
found a mean spectral index (7) = —0.09, with a dis- 
persion a = 0.28. Our results therefore confirm the 
steepening of extragalactic radio source sp ectra above 
~ 70 GHz, suggested by various d ata sets (jSadler et al.l 
20081; iGonzalez^Nueyo et all |200Sl : iMarriage et all l201ll ; 
Massardi et al. 2010l ) and then confirmed by analys es 
of different samples (|Planck Collaboration et al.l l2011bl lcl; 
iGiommi et al.l 1201 ll ) of bright extragalactic sources ex- 
tracte d by the Planck Early Release Com pact Source Cata- 
logue jPlanck Collaboration et alj|2011al ). 

This average steepening observed in the spectra of 
bright extragalactic s ources, which are mainly blazars in 
this frequency range (I Angel fc Stockman! [1980 1 . has been 
recently interpreted by iTucci et al.l (|201 lh in terms of the 
"break" frequency, foreseen in the emission spectra of 
blazars by classic physical models for the synchrotron emis- 
sion in the inner part of i nhomogenou s conic al jets. Remark- 
ably, this new model by ITucci et al.l (|201ll ) is able to give 
a good a fit not only to the source number counts here 
presented (see Figure [3] § 14. 1[) but also to almost all the 
published data on number counts and spectral index distri- 
butions of bright extragalctic radio sources up to, al least, 
200-300 GHz. 
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